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Abstract 

Nonstoichiometric perovskite-related oxides (such as ferrites and cobaltites, etc.) 
are characterized by fast oxygen transport at ambient temperatures, which relates 
to the microstructural texturing of these materials, consisting wholly of nanoscale 
microdomains. 

We have developed an inhomogeneous diffusion model to describe the kinetics 
of oxygen incorporation into nanostructured oxides. Nanodomain boundaries are 
assumed to be the high diffusivity paths for oxygen transport whereas diffusion 
into the domains proceeds much slower. Using Laplace transform methods, an exact 
solution is found for a ramped stepwise potential, allowing fitting of the experimental 
data to theoretical curves (in Laplace transforms) . A further model generalization is 
considered by introducing additional parameters for the size distribution of domains 
and particles. 

The model has been applied for qualitative evaluation of oxygen diffusion param- 
eters from the data on wet electrochemical oxidation of nano-structured perovskite 
SrCoo.5Feo.2Tao.303_y samples. 
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1 Introduction 



At high temperature grossly nonstoichiometric perovskite related oxides are 
solid state solutions with a high concentration of structural defects (oxygen 
vacancies, dopant ions, etc.). As temperature decreases, structural ordering 
occurs, producing either superstructures, in which defects are assimilated as 
structural elements, or complex microstructures with a high density of ex- 
tended defects [1]. The formation of the latter microstructures may be de- 
scribed using the concept of oversaturated solid solutions having a miscibility 
gap: oversaturation due to cooling causes unmixing of the initially homoge- 
neous solid solution. The unmixing can proceed either by spinodal decomposi- 
tion or homogeneous/heterogeneous nucleation [2]. Stabilization of nano-sized 
domains occurs due to a decrease in mobility of the domain boundaries at 
T < 0.3 T m . The mechanisms that are responsible for the enhanced thermal 
stability against microdomain growth are: (a) solute drag, (b) grain bound- 
ary segregation, (c) pinning by secondary phases and (d) chemical ordering 
[3]. Another approach to the explanation of intrinsic inhomogeneity of non- 
stoichiometric and doped oxides related to a strong tendency toward phase 
separation in systems with strongly correlated electrons. Electron conducting 
oxides with transition metals having different charges (for example, Mn 3 / 4+ , 
Co 3 / 4+ , Cu 2 / 3+ , etc) and magnetic states can break into a stable state made out 
of nanoscale coexisting clusters because of complicated electron-phonon and 
magnetic interactions [4] . The Gibbs energy of such a submicro- heterogeneous 
(nanostructured) system is essentially influenced by the nature of the interface 
between matrix and nucleus. This interface may be coherent, semi coherent 
and incoherent. If the interface is coherent or semicoherent the system behaves 
as a homogeneous single phase in terms of diffraction. Microdomain-textured 
oxides exemplify such systems rather well [1] . As shown in [5] , unmixing of non- 
stoichiometric and doped oxides with perovskite-related structures produces 
microdomains of 5-50 nm in size. 

Nanostructured materials may show enhanced ionic conductivity owing to the 
high density of interfaces which are enriched by defects due to the formation 
of space charge layers [6] or reduction of coordination number of cations in 
the vicinity of interface [7]. 

Therefore, nonstoichiometric doped oxides with microdomain texture may pro- 
vide conductive channels for enhanced oxygen diffusion along the microdomain 
boundaries [8,9]. Materials with high oxygen mobility are of interest for ap- 
plications such as oxygen electrolytes, electrodes, sensors and membranes for 
oxygen separation [8]. 

Numerous experimental proofs of enhanced diffusion along grain boundaries 
(including radiography, [10,11,12]) have initiated theoretical investigations 
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of inhomogeneous diffusion in polycrystalline materials [13,14,15]. The first 
mathematical description of fast diffusion along grain boundaries accompa- 
nied by penetration into the grains by means of bulk diffusion (without phase 
transition though) is shown in [13,16]. In [17,18], a mathematical model of the 
two-phase oxidation of MDT nonstoichiometric perovskite has been consid- 
ered. The crystal is conceived as an assembly of parallel domain boundaries, 
which are considered to be high diffusivity paths. Oxidation proceeds as a 
result of rapid oxygen transport along the domain boundaries and two-phase 
oxidation of the microdomains. 

Unlike the previous papers, here we assume that oxidation of microdomains 
occurs as a result of one phase reaction. We employ a model for oxygen trans- 
port in the MDT perovskites in the course of single phase oxidation, assuming 
that diffusion along the nanodomain boundaries is much faster than diffu- 
sion into the domains. Since nanostructured perovskites exhibit fast oxygen 
transport even at room temperature, the model was applied for the evalua- 
tion of oxygen diffusion parameters from data of wet electrochemical oxidation 
of nano-structured perovskite SrCoo.5Feo.2Tao.30 3 _ y samples in 1M KOH at 
ambient temperatures. 

Potentiostatic method [19,20] is considered to be a convenient tool for studying 
oxygen transport in perovskites (where electron conductivity far exceeds the 
ion one). In the present study we have used a version of this method. 



2 Experimental 

Kinetic studies were carried out for SrCoo.5Feo.2Tao.303_ y nonstoichiometric 
perovskite. The samples were synthesized by solid state reaction from the 
corresponding metal oxides and carbonates with preliminary homogenization 
of the starting materials in a planetary ball mill. A stoichiometric mixture 
of the powders was calcined at 900 °C, pressed in pellets, annealed in air at 
1400 °C for 6 hours and cooled in the furnace. 

For kinetic studies, the sintered sample was annealed at 950 °C in a quartz 
ampoule in dynamic vacuum (P ~ 10 3 Pa). After that, the quartz ampoule 
was closed and placed in water for rapid quenching of the sample to room 
temperature. 

Experiments were performed at room temperature (25 °C) in potentiostatic 
mode (three electrode cell, 1M KOH electrolyte, Hg/HgO reference electrode) 
with working electrodes of polycrystalline material (17-18 mg) pressed into 
Pt grids along with 1 wt.% of Teflon and 15-20 wt.% of acetylene black. The 
working electrode was placed in a cell at the set temperature and was main- 
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Fig. 1. Applied potential (a) and current (b) versus time. 

tained up to an equilibrium potential Uq. At t — a voltage pulse between the 
working electrode and reference electrode has been applied. Re-equilibration 
through diffusion of oxygen into the working electrode material takes place. 
The rate of re-equilibration can be measured easily by monitoring the electric 
current passing through the cell, see Fig. 1. 



3 Theoretical model 



Similar to the model of spheres filling a half-space (grains in a thick metal 
plate, see [14]), we consider a model, where small spheres are inside large ones 
(nano-domains inside the perovskite powder particles). The simple geometry 
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of our model allows us to derive an exact solution for the current passing 
through the cell using the Laplace transform method. 



Inverse Laplace transformation to determine the original function is not possi- 
ble in analytical form (in analytical functions; this is the reason why one usu- 
ally analyzes only the asymptotic current behavior at short and long times). 
Though there exist computational methods to do reverse Laplace transfor- 
mations [21], they require intensive calculations (double Fourier transforms). 
Therefore, we have chosen another way of model comparison with experiment. 
After we perform numerical Laplace transformation of measured current func- 
tion, Jexp(i), (as well as potential function, if it is measured and is not given 
"analytically" with potentiostate) we adjust transform, J exp (p), with model 
curves J t h(p) (which depend on the applied potential profile). 

However, computer inversion of Laplace transform may be useful for verifying 
the best approximation, or for the better choosing among several minima, if 
they occur. 

A two-level model of spheres is used to describe the current passing through 
an electrochemical cell under potential of a chosen profile. 



3. 1 Diffusion in a spherical particle 

The diffusion equation for a spherical particle, when the space derivatives 
(Laplacian) are reduced to differentiation over the radius, is written as follows: 



Here c(t, x) is the concentration of the diffusing species, R the particle radius, 
and x the radial coordinate, < x < R. For initial and boundary conditions 
we have 

c(t < 0, x) = , c(t, x = R) = c R (t) . 
Using the Laplace transform of c(t) 



d(x c) 
dt 



= D 



d 2 



2 (x c) 
dx 2 



(1) 



oc 
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turns partial derivatives (1) into an ordinary differential equation 

d 2 (xc) p 



dx 2 " D 



xc. (2) 



Applying boundary condition, c(p,R) = c R (p), one can easily find the exact 
solution of equation (2), 

x - / x kR sinh(/tx) / — — 

c(p,x) = c R (p) K = ^pD. 3 

Sinn (/cit) kx v 



Diffusion current (and electric one, if diffusing component carriers charge q), 
passing through the whole particle boundary, may be expressed as follows: 

R dc(t, x) 



J(t)=4nqJ^lx 2 dx, 



J(p) — Anq ^ R f ^' f smh(Kx) x dx . (4) 



o 



One may expect that the diffusion coefficient is independent of concentration, 
provided that the oxygen content in the perovskite lattice changes only slightly, 
e.g., Ay < 0.1. For this purpose the amplitude of potential pulse, Ua, applied 
to the cell, should be small enough (this is determined by the ratio of sample 
mass and charge). We assume that the access of oxygen ions to the particles 
is not limited and, hence, there is a linear dependence between concentration 
and the potential 

c fl (f) oc U(t) or c R (t)/c A = U(t)/U A . 



In experiments we used potential steps of the following form (see Fig. la): 



U(t)/U A = H(t;h,t ) = { 



0, t < 

h + (l- h)t/t , 0<t<t ; 

1, t > t 



this gives the Laplace transform (0 < h < 1): 

pc R (p)/c A = h + (l-h)(l- e-*>)/(pt ) ( ^ 1) • (5) 
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Fig. 2. Regions of slow (white) and fast (grey) diffusion 



Introducing value Q = J(0) = 47ri? 3 gc^/3 for a total charge passed through 
the cell, and using equations (3), (4) and (5) we obtain (r = R 2 /D is a char- 
acteristic diffusion time for a sphere of radius R; ^Jpf = kR) 

J(p)IQ=pc R (p)lc A F(^¥), (6) 
where F(y) = - 1 ( E=S 1 - + -..). (7) 



3.2 Inhomogeneous diffusion model: spheres of two types 

Let us imagine that particles in the sample are non-uniform (see Fig. 2) and 
contain domains of a typical size r, where diffusion rate is small, D ± , and well 
as regions neighboring with domains boundaries, where diffusion rate is high, 
D 2 . More exactly, let us consider D 2 as an effective diffusion coefficient for a 
porous particle, obtained, when all domains are removed, as if they are non- 
permeable for diffusion, i.e. Di = 0. This effective coefficient may change not 
only with temperature, but also with parameters related to the nano-domain 
structure, such as oxygen content, quenching rate, etc. 

If diffusion coefficient Di is small but not equal to zero, then diffusion equa- 
tion should reflect oxygen incorporation inside the domains. One may assume 
that the concentration along the boundary of each domain is approximately 
the same (since transport to the boundary is fast) and is equal (or at least 
proportional) to the current local concentration: 

ci r (t) = c 2 (x,t) . 

Therefore, we may write effective diffusion equation in a heterogeneous particle 
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instead of equation (1): 

-d- D2 i^-) = - ndJi ^ x)/q - 



Here rid is the concentration of domains, Ji/q the diffusion current into a 
domain, v\ = nd47rr 3 /3, v 2 — 1 — Vi are volume fractions for regions of 
slow and fast diffusion, respectively. After Laplace transformation we obtain 
(taking into account equations (5), (6)): 

^? 1 = ^ 2 (l + aF( v ^T)), (8) 



where a — v\jv 2 is the ratio of above mentioned volume fractions (see Fig. 2). 

As it is mentioned in [14], equilibrium concentrations (of diffusant) along the 
boundaries and inside the domains (i.e. grains) may not coincide. Therefore, 
dimensionless parameter a may include not only geometry but also the con- 
centration factor. 

One may derive the final expression for the current using equation (8) (it is 
necessary to remember that the particle is inhomogeneous) : 



Q 



x=R 

PC2R(P) 1 + ocf{pr 1 ) 



CA 



1 + a 



/(pr 2 [l + a/(p7i)]). 



(9) 



Here notations f(y) = F( v /y), T\ = r 2 / Di, t 2 = R 2 / D 2 (see equations (5)-(7)) 
are introduced. 



3.3 Assembly of particles, particle size distribution 



If powder particles are of different sizes, then the current expression (6) and 
expression (9) should be integrated with the function of volume (mass) distri- 
bution of particles, M(R) (norm per unit): 



J(p)/Q = ^M [ F( K R)M(R)dR. 
c A J 
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If the 'relative' dispersion of the distribution, 
A R = ((R-R ) 2 )/R 2 (R =(R)), 

is small, we leave the first correction term, changing functions in (9) or (6) as 
follows (prime means differentiation by argument): 

F(kR) -> F*(kR ) = F(kRo) + (kR ) 2 F"(kR )A r /2 , (10) 

where y 2 F"(y) = 6(1 + y 2 / sinh 2 y)(l + yj tanhy) - 24/y 2 . 

One may obtain this correction by e.g. integration with formal distribution 
(delta functions): 

M(R) = S(R - Ro) + 5"{R - R ) A R R 2 /2 , 

which has proper first two moments. Saddle-point integration with the Gaus- 
sian distribution M(R) oc exp(— (R/Ro — 1) 2 /2A R ) gives the same result. 

In the similar manner one may correct for the size dispersion of nano-domains, 



4 Experimental results and computation 

According to the X-ray diffraction (Fig. 3) and iodometric titration data, 
quenched samples of SrCoo.5Feo.2Tao.30"3_ y had the expanded cubic unit cell 
parameters (from 3.934 A to 3.952 A) and reduced oxygen content (from 2.92 
to 2.70) with respect to slowly cooled ones. 

XRD suggests single phase behaviour, while the HREM data shows that the 
samples possess a micro-domain texture (Fig. 4) with a typical domain size of 
about 10 nm. 

The mechanism of electrochemical oxidation of SrCoo.5Feo.2Tao.30"3__ y has been 
studied with chronopotentiometry combined with in situ X-ray diffraction. 
The results of evaluation of the in situ data are shown in Fig. 5. Monotonous 
changes in the unit cell parameters and potential vs. inserted oxygen (x) (x = 
n/2, where n is a charge transfer) provide evidence of one phase mechanism 
of oxidation [22]. 

For kinetic studies, the samples were powdered, and particle size distribution 
was measured (see Fig. 6). 
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Fig. 3. Diffraction patterns of samples SrCoo.5Feo.2Tao.303_ y : (1) as sintered, (2) 
after annealing at 950 °C and quenching in vacuum (P ~ 10 3 Pa). 




Fig. 4. High resolution image of SrCoo.5Feo.2Tao.3O2. 7 sample possessing microdo- 
main texture. 

Fig. 1 shows U(t) and J cxp (t) functions versus time, obtained in one of ex- 
periments. The time of current measurement averaging was 1 second. This 
obstacle somewhat distorted the current plot at short times. 

Computations were performed in MatLab. For the data, presented in Fig. 1, 
potential step parameters are the following (see equations (5) and (9)): 



h = 7/55.25; t_0 = 958; %% sec 
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Fig. 5. In situ anodic oxidation of SrCoo.5Feo.2Tao.302.7+ x : (a) potential E vs. charge 
transfer n, (b) change of unit-cell parameters with n, n/2 = x. 
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Fig. 6. Mass distribution of particles over size. 

After J ex p(p)/Q was calculated (exponential slope at long times (t >3000 s) 
was considered analytically), correction for the time of micro ammeter inte- 
gration (1 s) was introduced: 

Jep = Jep.*p./(1 - exp(-p)); 

Then optimization over three parameters entering the model function of cur- 
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Fig. 7. Normalized Laplace transform of experimental data (solid; Q = 0.44 C) 
and fitting model curves for inhomogeneous (dot) and homogeneous (dash; a = 0) 
model. 

rent (equation (9)) was done: 
parO = [taul,tau2, alpha]; 

options = optimset('TolFun',le-10,'TolX',le-6, ... 

'MaxFunEvals',le3,'MaxIter',le3); 
[parl,fval,exitflag,output] = ... 
fminsearch('diff_fun',parO, options, p,Jep,handle); 

function ff = diff_fun(param,p,Jep, handle) 

ff = sum(abs(Jtp - Jep)); 

Computing results are shown in Fig. 7, model parameters for theoretical curve 
(dotted line) being as follows: 

n = 49 s, t 2 = 19 • 10 3 s, a = 4.8 . 



Using values R = 32 /im, r ~ 10 nm, one may estimate diffusion coefficients 



D 2 = 5 • 10~ 10 cm 2 /s, Di = 2 ■ 10~ 13 cm 2 /s 
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Let us note that fitting by the homogeneous model gives essentially larger 
errors (see dashed curve in Fig. 7), and yields another value: r 2 = 15 • 10 4 s. 
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5 Discussion and conclusions 



As we have already noted, nonstoichiometric and/or doped oxides are high 
temperature systems, since solid solutions are stable at relatively high tem- 
peratures only At lower temperatures the high free energy related to the high 
concentration of defects (oxygen vacancies, guest ions) can be reduced either 
by defects ordering and their localization as structural elements (superstruc- 
tures formation), or by their elimination or precipitation (nanostructuring, 
microdomain texturing) [1]. The difference between the MDT phases and or- 
dered superstructures is that in MDT oxides part of disorder is already accu- 
mulated at low temperatures as excess interfacial energy (disordered lattice 
in the vicinity of interfaces). Therefore, the lower is working temperature, 
the more essential is the MDT oxide advantage with regard to oxygen trans- 
port. That is why MDT oxides are promising materials for low and moderate 
temperature oxidative catalysis, oxygen permeable membranes for partial ox- 
idation of hydrocarbons, SOFC's electrodes, sensors etc. 

With our developed model, we have shown that nanostructuring allows inho- 
mogeneous oxygen diffusion in domains and along the interfaces; the difference 
in diffusion coefficients may attain several orders of magnitude. The most im- 
portant issue is that stationary flux of oxygen ions through a MDT perovskite 
membranes at temperatures lower than the point of order-disorder transition 
is determined mostly by coefficient D 2 . 

Since activation energy E a , necessary for the oxygen ions migration along the 
grain boundaries (g.b.), may be essentially lower than that for the bulk (b.) 
diffusion, i? a (g.b.) ~ l/2E a (b.) [23], oxygen permeable membranes made of 
MDT oxides are able to provide several orders of magnitude higher oxygen 
fluxes at the working temperatures below the order- disorder transition point. 
Let us mention though that nanostructuring is not the sufficient condition for 
intensive oxygen transport along the interfaces. For example, alkali/rare earth 
manganites possess a well-developed microstructure, but oxygen mobility in 
them is very low, whereas alkali/rare earth ferrites, and in particular cobaltites 
show unusually high oxygen transport [8]. Apparently, the transport properties 
of oxides are influenced by the interfacial energy, which is dependent on a 
set of parameters (M-0 bond strength, coordination of cations, their charges 
and size, electronic configuration, magnetic state of cations, etc.). In case if 
interfacial energy is high enough (disordered interface), one may expect high 
permeability of oxygen. If interface is ordered (low interfacial energy), most 
likely it will be a barrier for the ion diffusion. 

We shall further study the kinetics of MDT perovskites oxidation at various 
temperatures to determine the activation energy for oxygen ions migration 
along the interfaces and inside the domains and develop the method for mea- 
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suring oxygen diffusion parameters in the nanostructured materials. Compar- 
ing the data obtained at low temperature by means of relatively simple wet 
electrochemical technique with the data related to oxygen permeability at high 
temperatures, we might define the mechanism of oxygen transport in MDT 
oxides, elucidate the factors determining the high values of oxygen fluxes in 
these materials and might open strategies to develop new oxygen-conducting 
materials operating at moderate temperature. 
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